Strongly Enhanced Vortex Pinning by Conformal Crystal Arrays 
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Conformal crystals are non-uniform structures created by a conformal transformation of regular 
two-dimensional lattices. We show that gradient-driven vortices interacting with a conformal pinning 
array exhibit substantially stronger pinning effects over a much larger range of field than found for 
random or periodic pinning arrangements. The pinning enhancement is partially due to matching 
of the critical flux gradient with the pinning gradient, but the preservation of the sixfold ordering in 
the conformally transformed hexagonal lattice plays a crucial role. Our results can be generalized 
to a wide class of gradient-driven interacting particle systems such as colloids on optical trap arrays. 
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One of the most important problems for applications 
of type-II superconductors is how to create high critical 
currents or strong vortex pinning over a wide range of ap- 
plied magnetic fields^. For over sixty years, it has been 
understood that the ground state vortex structure is a 
hexagonal lattice 2 -, so many methods have been devel- 
oped to increase the critical current using uniform pin- 
ning arrays that incorporate periodicity to match the vor- 
tex structure^—. The pinning is enhanced at commensu- 
rate fields when the number of vortices equals an integer 
multiple of the number of pinning sites, but away from 
these specific matching fields, the enhancement of the 
critical current is losti^. Efforts to enhance the pinning at 
incommensurate fields have included the use of quasicrys- 
talline substrates^ or diluted periodic arrays^—, where 
studies show that new types of non-integer commensu- 
rate states can arise in addition to the integer matching 
configurations. Hyperbolic tessellation arrays were also 
recently considered 21 . 

Part of the problem is the fact that under an applied 
current, the vortex structure does not remain uniform 
but instead develops a Bean-like flux gradient^ 2 .: the vor- 
tex density is highest at the edges of the sample when 
the magnetic field is increased, and highest in the center 
of the sample when the magnetic field is removed and 
only trapped flux remains inside the sample. As a con- 
sequence, uniform pinning arrays generally have a por- 
tion of the pinning sites that are not fully occupied, sug- 
gesting that a more optimal pinning arrangement should 
include some type of density gradient to match the crit- 
ical flux gradient. Here we show that a novel type of 
pinning structure, created using a conformal transfor- 
mation of a uniform hexagonal lattice, produces a much 
higher critical current over a much wider range of mag- 
netic fields than any pinning geometry considered up un- 
til now. Conformal crystals not only have a density gra- 
dient, but also preserve aspects of the hexagonal order- 
ing naturally adopted by the vortex lattice. The pinning 
enhancement we find is substantial and will be very im- 
portant for a wide range of superconductor applications 
and flux control. Our results can also be important for 
stabilizing novel self-assembled structures created using 
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Figure 1: A conformal transformation is applied to the semi- 
annular section of a regular hexagonal lattice shown in (a) to 
create the conformal crystal structure shown in (b)^. Points 
A — F in (a) are mapped by the transformation to points 
a — f in (b) respectively. The straight contour lines connect- 
ing nearest neighbor lattice points in (a) are bent into arcs 
in (b), but the local six- fold ordering of the lattice points is 
maintained. Pinning sites are placed at the vertices formed 
by the intersections of the contour lines in (b). 



density gradients, such as in colloidal systems 2 ^. 

To examine vortex pinning by conformal crystalline 
substrates, we numerically compute the magnetization 
M, which is a measure of the pinning strength. We com- 
pare the effectiveness of the conformal pinning to that 
of random and periodic pinning arrays with the same 
number of pinning sites. The conformal array produces 
significantly stronger pinning over a much larger range of 
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fields than these other arrays, with the exception that the 
pinning by periodic arrays is strongest only in the vicin- 
ity of the matching fields 14 . By comparing our results 
to those obtained for random pinning with an equivalent 
gradient, we demonstrate that it is not simply the non- 
uniform density of the conformal array that gives rise to 
the pinning enhancement; rather, it is the preservation 
of the hexagonal ordering of the pinning by the confor- 
mal map that plays the most critical role. Another fea- 
ture that gives the conformal pinning an advantage over 
periodic pinning for incommensurate fields is that even 
though local hexagonal ordering is present, the overall 
arc-like arrangement of the pinning sites (apparent in 
Figs.QJb) and[2|) prevents the formation of straight-line 
channels along which vortices can easily flow. In addi- 
tion to vortex pinning in type-II superconductors, the 
effects of conformal crystalline substrates on ordering or 
dynamics of a monolayer of particles could also be stud- 
ied for vortices in Bose-Einstein condensates on optical 
lattices 2 ^ or colloidal particles on optically created sub- 
strate arrays^ 3 -. The enhanced pinning also suggests that 
conformal arrays could be used to increase friction for 
particle-surface interactions. 

Conformal crystals are a class of two-dimensional (2D) 
structures created by the application of a conformal 
(angle-preserving) transformation to a regular lattice in 
the complex pian o 25 ' 26 . Figure Q] illustrates a conformal 
crystal obtained via the transformation of a hexagonal 
lattice. The contour lines connecting nearest neighbors, 
which are straight lines for the original hexagonal lat- 
tice, are bent into arcs but still cross at angles of 7r/3, 
preserving the sixfold coordination of individual pinning 
sites in spite of the clear density gradient. To create 
a pinning lattice, we place pinning sites at the vertex 
locations where the contour lines intersect. Conformal 
crystal structures have been studied experimentally for 
repulsively interacting magnetic spheres confined to a 2D 
container that is tilted so that the gravitational force on 
the particles produces a mechanically stable but nonuni- 
form crystal 2 ^. Other systems where conformal crystals 
arise include foams under an external fiel d 28 ' 29 and large 
arrays of classical Coulomb charges in confined circular 
potentials where local triangular ordering occurs along 
curved lattice lines 3 -^. 

Simulation- We conduct a flux gradient density simu- 
lation of the type previously used to study vortex critical 
states and magnetization with rando m 31 ' 32 and periodic 
pinning arrays 1 ^. Figure [5] shows the outer pin-free re- 
gion surrounding a central pinned region that consists 
of two conformal crystals placed with their highest den- 
sity regions adjacent to the pin-free region. Details on 
the construction of the conformal crystal are given in the 
supplemental material^!. We use periodic boundary con- 
ditions in the x and y-directions and consider a 36A x 36A 
system with pinned region extending from x — 6 A to 30 A, 
where A is the penetration depth. This geometry was pre- 
viously shown to be large enough to capture accurately 
the behavior of the magnetization curve s 14 ' 31 ' 32 . 
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Figure 2: The gradient-driven sample geometry consists of 
two conformal crystals facing each other in the pinned region. 
Open circles are pinning site locations. Vortices are added to 
or subtracted from the pin-free region labeled "A" on the left 
and right sides of the image. 



The dynamics of vortex i are obtained by inte- 
grating the overdamped equation r](dRi/dt) — F™ + 



?] is the damping constant which is set equal to 
The vortex-vortex interaction force is F'" v = 
where K\ is the modified 



unity. 

Yh=i s l s j F K 1 (R ij /X)'R l: j, 



Bessel function, R, is the location of vortex i, JL» = 



R, — R-; , R, 



(R 



-j i' - ~ Rj)/ R iv F o = ^(FMoA 3 , and tp 
is the elementary flux quantum. The sign prefactor is 
+1 for a vortex and —1 for an antivortex. The pinning 
sites are modeled as N p non-overlapping parabolic traps 



with 17 = Tti(F p R p jR P )Q((Rv-R%)im^ where 
i?| is the location of pinning site k, R p ik = |R, — R?|, 

Rf fe = (Rj — TlT)/R? k , is the Heaviside step function, 
R p is the pinning radius that we fix to R p = 0.12A, and 
F p is an adjustable parameter controlling the strength 
of the pinning force. All forces are measured in units of 
Fo and all lengths in units of A. The external field is 
measured in units of H ( j >1 the field at which the average 
unit density of vortices equals the average unit density of 
pinning sites. 

To perform a complete field sweep, we begin with 
zero vortex density and then quasistatically add vortices 
in the unpinned region (labeled "A" in Fig. [5]) at ran- 
domly chosen nonoverlapping positions until the desired 
maximum external field value is reached. As the vor- 
tex or antivortex density builds up in the pin-free re- 
gion, the vortices or antivortices drive themselves into 
the pinned region due to their own repulsive interactions, 
creating a gradient in the flux density within the pinned 
regio n 14 ' 31 ' 32 . We then remove vortices from the pin- 
free region until the vortex density in this region reaches 
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Figure 3: The magnetization M vs H/H^, where is the 
matching field at which the vortex density equals the average 
pinning density, in samples with n v = 1.0 and F v = 0.55. 
Outer dark (black) curve: a sample with a conformal pinning 
array (CPA). Inner light (red) curve: a sample with a uni- 
formly dense random arrangement of pinning sites. The CPA 
produces a much higher value of M at all fields. 



zero. To reverse the field, we add antivortices, which 
repel each other but are attracted to vortices, to the 
unpinned region. When a vortex and antivortex come 
within a small distance (0.3A) of each other, they are 
both removed from the system to simulate an annihila- 
tion event. To complete an entire magnetization loop, 
we continue to add anti-vortices until the external field 
reaches its most negative value, and then remove antivor- 
tices from the pin-free region to bring the external field 
back up to zero. The average magnetization M is calcu- 
lated as the difference between the flux density H in the 
unpinned region and the density B in the pinned region, 
M = -(1/4ttV) f(H - B)dV, where V is the sample 
area. The critical current J c is proportional to the mag- 
netization. 

Results- In Fig. [3] we plot the complete hysteresis loop 
M vs H/H^ for the conformal pinning array (CPA) and 
the uniformly dense random pinning array. Each sample 
contains the same number N p of pinning sites with an av- 
erage pinning density of n p = 1.0 and with F p = 0.55. We 
find that M is much higher at all fields for the CPA than 
for the random pinning, and that at its highest point, 
M for the CPA is almost four times higher than for the 
random pinning. Although the CPA has local triangular 
ordering, we observe no peaks or other anomalies in M 
at integer matching or fractional matching multiples of 
H/H<f, of the type found for uniformly dense periodic pin- 
ning arraysii. The flux profiles plotted in Suppl. Fig. 1— 
show that the random array produces a Bean-like profile 
that becomes shallower as H increases. In contrast, at 
higher fields the CPA does not have a uniform flux gra- 
dient but instead develops a double slope profile, with 
a larger flux gradient near the edge of the sample and 
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Figure 4: (a-c) M vs H/H^ plots for samples with F p = 0.55 
and average pin density n v — 1.0. (a) Dark (black) curve: the 
CPA system from Fig. |3j light (red) curve: a uniform square 
pinning array, (b) Dark (black) curve: the CPA system; light 
(red) curve: a uniform triangular pinning array. M is higher 
for both the uniform periodic pinning systems than for the 
CPA at the matching field near H/Hj, = 1.0, but the CPA 
has the highest value of M at all other fields, (c) Outer dark 
(black) curve: the CPA system; middle light (green) curve: a 
random arrangement of pinning sites with an equivalent pin- 
ning gradient; inner light (red) curve: the uniform random 
array. Here the gradient in the random pinning increases M 
compared to the uniform random pinning by only a small 
amount for most fields, (d) Mhw, the half- width of the mag- 
netization loop at H /H^ = 1.0, for the CPA (circles), random 
array with CPA-equivalent pinning gradient (triangles), and 
uniform random pinning array (squares) for varied F p . There 
is a consistent enhancement of M in the CPA. 



a much shallower or nearly flat flux profile in the cen- 
ter of the sample. As H increases, the sharper slope 
region decreases in width and is replaced by the shallow 
slope region. We have examined a range of pinning forces 
and find that the magnetization for the CPA is consis- 
tently higher than for random pinning arrays, as shown 
in Fig. Hid) where we plot Mhw, the half- width of the 
magnetization loop at H/H<p = 1.0, at varied F p for these 
array types. 

We next address whether the conformal pinning arrays 
produce higher pinning compared to other previously 
studied types of non-random pinning arrays. In Fig.[3ta) 
we plot M vs H/H^ for the CPA and a square pinning 
array with the same pinning density and strength, and 
in Fig.0Jb) we plot the same quantity for the CPA and a 
triangular pinning array. In both cases, M for the CPA 
is higher over most of the range of H/H^ except at the 
first matching field, where the periodic pinning arrays 
produce a higher value of M due to a commensurability 
effect. This shows that although a periodic pin structure 
can strongly enhance the pinning, this enhancement oc- 
curs only for a very specific matching field. In contrast, 
the CPA produces a significant enhancement of the pin- 
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Figure 5: (a) The high field behavior of M vs H/H^, for 
samples with F p = 0.55 and n p — 1.0 on the initial ramp 
up only. Dark (black) curve: CPA; light (red) curve: uniform 
random array. Although M is higher for the CPA than for the 
random array over the entire range of fields, there is a drop in 
M for H/H<j, > 2.0 in the CPA. (b) The fraction of occupied 
pinning sites P vs H/H^, for the CPA (dark black curve) and 
uniform random pinning (light red curve). P is initially higher 
for the CPA but drops at H/H,f, = 2, coinciding with the drop 
in M in (a), while P for the random pinning sample increases 
monotonically with increasing field. 

ning over a very broad range of fields, extending well 
above the first matching field. 

Since the CPA has a pinning gradient, it could be pos- 
sible that any type of pinning array with an equivalent 
gradient would also exhibit a pinning enhancement com- 
pared to uniform pinning arrays and could be just as effi- 
cient at pinning as the CPA. We hnd that this is not the 
case. In Fig. 0{c) we plot M versus H/H^ for the CPA 
and a random pinning array with an equivalent pinning 
gradient and the same average pinning density. For com- 
parison we also show the uniform random pinning array. 
The random pinning with a gradient exhibits a modest 
enhancement of M compared to the uniform random pin- 
ning array; however, both arrays give M values that are 
significantly smaller than the CPA for all but the very 
lowest fields. This result indicates that it is the other 
properties of the CPA such as the preservation of local 
sixfold ordering, and not merely the pinning gradient, 
that are responsible for the enhanced pinning. The arch- 
ing structure of the CPA suppresses certain modes of 
vortex motion. For example, for random pinning arrays 
it has been found that certain regions where the pinning 
density is slightly lower lead to the formation of persis- 
tent river-like flow patterns 3 ^; these structures are unable 
to form in the CPA. For periodic pinning arrangements, 
easy vortex flow occurs along the symmetry directions 
of the pinning arra y 33 i 34 , causing a drop in the critical 
current above a commensurate field. Such motion is sup- 
pressed in the CPA due to the arching structure. 

In Fig. EJa) we compare M versus H/ for the CPA 
and uniform random arrays in the first quarter of the 
magnetization loop over a much larger range of fields up 
to H/Htf, = 4.0, and in Fig.jSfb) we plot the correspond- 



ing pin occupancy P, which is the fraction of pinning 
sites occupied by vortices. The enhanced pinning for the 
CPA is the most pronounced below H/H^ — 2.0. The en- 
hancement decreases above this field but remains larger 
than the random pinning array for all fields. For the ran- 
dom pinning array, P monotonically increases over the 
entire range of H/H^. In contrast, after running well 
above the P value for the random pinning array at lower 
fields, P for the CPA rolls over and begins to decrease 
with increasing field above H/H<p w 2, correlated with 
the decrease in M . This is also the field at which the 
higher gradient region seen in Suppl. Fig. 1(a) disappears 
from the sampled. Just below this field, all of the pin- 
ning sites near the edge of the sample are occupied, and 
since these pinning sites are the most closely spaced, as 
additional vortices enter the sample, the vortex- vortex in- 
teractions become strong enough to push some of the vor- 
tices out of the pins, leading to the drop in P and M. For 
the random pinning array, there are always some empty 
pinning sites near the edge of the sample in places where 
two pins happen to be in close proximity, so that the 
vortex-vortex interaction energy would be prohibitively 
high if both pins were occupied simultaneously. As the 
field increases, these pinning sites gradually become oc- 
cupied. Even though P for the CPA falls below P for the 
random array at higher fields, the pinning enhancement 
remains significantly stronger for the CPA. 

In summary, we demonstrate strongly enhanced vortex 
pinning by a conformal crystal array of pinning sites. The 
conformal crystal is constructed by a conformal transfor- 
mation of a hexagonal lattice, producing a nonuniform 
structure with a gradient where the local sixfold coor- 
dination of the pinning sites is preserved, and with an 
arching effect. The conformal pinning arrays produce 
significantly enhanced pinning over a much wider range 
of field than that found for other pinning geometries with 
an equivalent number of pinning sites, such as random, 
square, and triangular. We show that the pinning en- 
hancement is not simply due to the pin density gradient, 
but is also due to the preservation of the sixfold coor- 
dination of the pinning sites and to the arching effects 
that prevent the formation of easy channels of vortex 
flow. The conformal pinning arrays do not show any 
pronounced enhancement at the matching fields of the 
type found for periodic pinning arrays; however, there is 
a drop in the magnetization at higher fields that is corre- 
lated with a drop in the occupancy of the pinning sites. 
Our results will be important for applications where high 
critical currents are required, and can be generalized to 
stabilizing or pinning other systems, such as colloidal par- 
ticles, that do not form uniform crystalline structures. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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SUPPLEMENTARY MATERIAL 

for "Strongly Enhanced Vortex Pinning by 
Conformal Crystal Arrays" 



of the main text, specified by 

< u < 7rr out and < v < r out In 



\ ' in 



A. Explicit Construction 



I. CONSTRUCTION OF THE CONFORMAL 
CRYSTAL 

To obtain the conformal crystal structure, we first sit- 
uate a regular hexagonal lattice in the complex plane: 
that is, we encode the location (x, y) of each lattice site 
into a complex number z = x + iy. With lattice constant 
6, the lattice sites will be located at 



z = m • (1 • b) + n 2 ■ (e ra/3 • &) 



(1) 



where n\ and n 2 are arbitrary integers. 

We then consider a semiannular region of this lattice 
as shown in Fig. 1(a) of the main text, defined by Imz > 
and n n < \z\ < r out , and apply a conformal (angle- 
preserving) transformation mapping z to w = u + iv, 
given by 



w = 1 In (laz) 

2a lot 



(2) 



where a = ^— . The semiannular region in the z-plane 
shown in Fig. 1(a) of the main text will be mapped to 
the rectangular region in the w-plane shown in Fig. 1(b) 




Figure 1: Data from the samples in Fig. 3 of the main text 
with n p = 1.0 and F p — 0.55. (a) The average flux density 
p v across the sample plotted for different times during the 
ramp up portion of the hysteresis curve for the conformal 
pinning array (CPA), (b) p v for the random pinning array. 
The random array exhibits a Bean profile while the CPA has 
two regions inside the sample: a higher flux gradient at the 
edges of the pinned region and a small gradient in the middle. 



The structure of the final conformal crystal is com- 
pletely determined by specifying the three constants 6, 
r out , and r- ln . The locations of the sites of the conformal 
crystal can be generated directly by the expressions 

'2ni + n 2 



U = r out 



{5- 



tan 



n 2 



V3 



v = r out In 



Tout 



bsfr 



n x n 2 



(3a) 
(3b) 



where (ni,n 2 ) range over all pairs of integers satisfying 
the two constraints, n 2 > and 



<b 2 



+ nin 2 ) < r 



2 

out ' 



To obtain a conformal crystal of length u max and width 
u max , one chooses r out and r- m to have the values 



Tout = 



u a 



^max 

r in = exp 



b can be chosen to obtain a desired final density p of 
lattice sites: 

, 2 1 - (n„Aout) 2 



pV3 ■ In (r ut/r in ) 

Alternatively, it can be set equal to a desired minimum 
distance between pinning sites in the final structure. 
(Note that in Fig. 1 of the main text, the lengths of arc 
ABC and line abc are equal.) 

Figure 1 of the main text was obtained using r out = 
36/tt, r in = we-^ 3 , and b 2 = (1 - exp (-2tt/3)) • 
(V3/7t). These parameter values were calculated from 
u m a X = 36, w max = 12, and p = 1. 



II. FLUX DENSITY PROFILES 

In Supplemental Fig. 1(a) we plot the flux density pro- 
files obtained by averaging the flux density in the y di- 
rection for the first ramp up of the field in the hysteresis 
loop for the conformal pinning array (CPA) illustrated 
in Fig. 3 of the main text, and in Supplemental Fig. 1(b) 
we show the flux density profiles for the uniform ran- 
dom pinning array from Fig. 3 of the main text. The 
random array produces a Bean-like profile that becomes 
shallower as H increases. In contrast, at higher fields the 
CPA does not have a uniform flux gradient but instead 
develops a double slope profile, with a larger flux gradi- 
ent near the edge of the sample and a much shallower or 
nearly flat flux profile in the center of the sample. As 
H increases, the sharper slope region decreases in width 
and is replaced by the shallow slope region. 



